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Abstract 

We investigate the QCD magnetic susceptibility x the finite quark-chemical potential ^ 0) 
and at zero temperature (T = 0) to explore the pattern of the magnetic phase transition of the 
QCD vacuum. For this purpose, we employ the nonlocal chiral quark model derived from the 
instanton vacuum in the presence of the chemical potential in the chiral limit. Focusing on the 
Nambu-Goldstone phase, we find that the magnetic susceptibility remains almost stable to ~ 200 
MeV, and falls down drastically until the the quark-chemical potential reaches the critical point 
fic ~ 320 MeV. Then, the strength of the x is reduced to be about a half of that at = 0, and 
the first-order magnetic phase transition takes place, corresponding to the chiral restoration. From 
these observations, we conclude that the response of the QCD vacuum becomes weak and unstable 
to the external electromagnetic field near the critical point, in comparison to that for vacuum. It 
is also shown that the breakdown of Lorentz invariance for the magnetic susceptibility, caused by 
the finite chemical potential, turns out to be small. 
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I. INTRODUCTION 



It is of great importance to understand nontrivial structures of the QCD vacuum, since 
it reflects strong nonperturbative fluctuations and lay the foundation for explaining low- 
energy phenomena of hadrons. In particular, the quark condensate characterises the non- 
perturbative features of the QCD vacuum and plays a role of an order parameter associated 
with spontaneous breakdown of chiral symmetry (SB^S) which is known to be essential in 
low-energy hadronic phenomena. The nature of the QCD vacuum can be examined under 
external fields such as electromagnetic (EM) field, since it reveals the EM properties of the 
vacuum. The magnetic susceptibility is defined as the vacuum expectation value (VEVs) in 
Euclidean space: 

= ieqF^^{iq^q)x, (1) 

where Cq denotes the electric charge of the quark and F^y represents the EM field-strength 
tensor. The {iq^q) stands for the quark (chiral) condensate taken as a normalization. It is 
natural to take it as a normalization for the magnetic susceptibility, since the VEV {iq^a^^q)F 
breaks chiral symmetry [T]. Recently, Braun et al. [2] have proposed that the x may be 
measured in the exclusive photoproduction of hard dijets + N qq + N. Theoretically, 
there have been already many works on it in the QCD sum rules and chiral effective models [3l 
in El El m [HI [9]. The magnetic susceptibility is also related to the photon distribution 
amplitude p]. 

While the magnetic susceptibility is relatively well understood in free space, it is not 
much known about its modification at finite density (finite quark-chemical potential, /i). 
For example, there are several interesting theoretical points as follows: First, in principle, 
the pattern of the magnetic phase transition of QCD vacuum can be revealed by studying 
the change of the magnetic susceptibility in medium [10]. Second, one can see how the 
vacuum reacts in the presence of the external EM source, especially, near the critical point. 
Third, the medium modification of the x ^lay play a critical role in explaining the exclusive 
photoproduction of the hard dijets in nuclei. Finally, this information may be important 
in understanding the EM features of neutron stars [H] . Moreover, since Lorentz invariance 
is broken in medium, we must consider separately the space and time components of the 
magnetic susceptibility at finite quark-chemical potential. In fact, Refs. [121 [131 E] have 
studied the importance of the breakdown of Lorentz invariance for the in-medium pion weak 
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decay constant F.„. In fact, Ref. [15] has calculated the magnetic susceptibility without 
the breakdown of Lorentz invariance taken into account. In addition, the chemical potential 
was treated perturbatively in Ref. (TS]. We will consider explicitly the breakdown of Lorentz 
invariance in the present work. 

In order to compute the x finite quark-chemical potential, we employ the framework 
of the yU-modified nonlocal chiral quark model (NLxQM) derived from the instanton vac- 
uum [16] . The NLxQM has several virtues: All relevant QCD symmetries are satisfied within 
the model. In particular, the instanton vacuum well realizes spontaneous breakdown of chiral 
symmetry via quark zero modes. Moreover, there are only two parameters: the average size 
of instantons (p ~ 1/3 fm) and average inter-instanton distance (-R ~ 1 fm), which can be 
determined by the internal constraint such as the self-consistent equation [TTl UHl [19] . These 
values for p and R have been supported in various LQCD simulations recently [20|, [211 122] . 
There is no further adjustable parameter in the model. Note that the inverse of the av- 
erage size of instantons, p^^ ~ 600 MeV, is regarded as the normalization point of the 
present model, so that we can calculate the dependence of the magnetic susceptibihty on 
the chemical potential to the critial point {p ~ 320 MeV) without any trouble. 

We concentrate in this work on the Nambu-Goldstone (NG) phase. In order to go beyond 
the NG phase, we need subtle techniques which we will explain briefly later. As a result, we 
find that the x remains stable up to p ^ 200 MeV, and falls down drastically to the critical 
quark-chemical potential, pc ~ 320 MeV and the strength of the x is reduced to be about a 
half of that at p = 0. The first-order magnetic phase transition takes place at this critical 
point corresponding to the phase of the chiral restoration. We will show that the response 
of the QCD vacuum becomes weak and unstable to the external EM source near the critical 
point. Moreover, the breakdown of Lorentz invariance turns out to be small in contrast to 
that in the case of the pion weak decay constant F^^ [12] . 

We organize the present work as follows: In Section II, we briefly explain the /x-modified 
NLxQM. In Section III, we compute the magnetic susceptibility x ^^id present the numer- 
ical results with discussions. The final Section summarises the present work and draws 
conclusions. The outlook of this work is also shortly discussed in this Section. 
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II. NONLOCAL CHIRAL QUARK MODEL AT FINITE QUARK-CHEMICAL 
POTENTIAL 



A modified Dirac equation with the finite quark-chemical potential n in (anti)instanton 
ensemble can be written in Euclidean space as follows: 

[tp-z^+An]^^;j' = K^^;j>, (2) 

where fia = (0,0,0,/i). Note that we consider here the chiral limit (m — > 0). The subscript 
(/)/ stands for the (anti)instanton contribution, and we use a singular-gauge instanton 
solution: 

^^^^^ = x^Jx^ + p^y 

where r/^'' and p denote the 'tHooft symbol and average instanton size (p ~ 1/3 fm), re- 
spectively. Since the instanton has its own spatial distribution with finite size (p) in prin- 
ciple, one needs to consider an instanton-distribution function when the integral over the 
collective coordinate for the instanton (instanton center, color orientation, and instanton 
size) is performed [171 [181 HH]- However, we assume here a simple 5-fucntion type for the 
insatanton-distribution function, i.e., 6{p — p)- In fact, the effect of the finite size is taken 
as 1/Nc corrections [HI [19]. We refer to Ref. [9| for the finite-size effect on the magnetic 
susceptibility in free space. The quark zero-mode solution then can be obtained by solving 
the following equation: 

[ip + A,j] = 0. (4) 

The explicit form of can be found in Refs. [HI [23]. Assuming that the low-energy 
hadronic properties are dominated by the zero mode, we can write the quark propagator in 
one instanton background approximately as follows: 

Siiix, y) = mx)ij{yy) = - E ""n^-^)^!! ^ ^ So{x, y) - ""n^^)^!! ^V) ^ (5) 

„ An + im tm 

where 5*0 is a free quark propagator, 5*0 = (if) — ifi)^^. Incorporating Eqs. ^ and ([s]), and 
averaging over all instanton collective coordinates, we can obtain the following expression 
for the quark propagator in momentum space: 



The momentum- and /i-dependent effective quark mass M{p, fi), corresponding to the Fourier 
transform of the quark zero-mode solution of Eq. (|4]), reads: 

M{p,fi) = Mo(/i) [{p + zfifij\p,fi)], (7) 

where Mq designates the constituent-quark mass as a function of fi, which will be determined 
self-consistently within the model. The analytical expression for ipip^fi) is also given in 
Ref. [16] . In actual calculations, we make use of a dipole-type parameterization, instead of 
using Eq. ([T]), for numerical simplicity: 



M{p,fi) ^ Moifi) 



2A^ 

{p + ifiy + 2A2 



(8) 



The cutoff mass A corresponds to the inverse of the average (anti)instanton size (~ 1/p), 
resulting in A ^ 600 MeV. We have verified that this simple parameterization reproduces 
Eq. ([T]) qualitatively well [21]. One can refer to Ref. [25j for more quantitative studies on 
the parameterization of M{p,fi). 

Now we are in a position to consider the effective partition function for Nf = 2 from 
the instanton vacuum configuration with the finite fi in the large Nc limit. Choosing the 
relevant terms for further discussion, it can be written as follows [16]: 



Zes = J dX Dip Dip^ exp J d'^xip^i^ - + im)tp + X{Y+ + ¥') - N\nX 



(9) 



where the flavour, colour, and Dirac spin indices are assumed tacitly. Note that this effective 
partition function has been constructed in such way that we can obtain the quark propagator 
given in Eq. ([s]). denotes the 2A^y-'tHooft interaction in the (anti) instanton background 
at finite /i. As for the parameters appearing in the second and third terms of Eq. (|9]), A 
plays a role of the Lagrangian multiplier, whereas indicates the sum of the number of 
the (anti)instantons. For more details on this partition function and its derivation, one may 
refer to Ref. [16], and references therein. 

We would like to explain briefly the phase structure obtained via the present model. 
Since both the quark- ant iquark and quark-quark interactions for Nf = 2 are attractive, 
two possible phases arise: The Nambu-Goldstone (NG) and colour-superconducting (CSC) 
phases. They are characterised by the QCD order paramters, i.e., the chiral and diquark 
condensates, respectively, or by the finite constituent quark mass Mq and diquark energy 
gap A, equivalently. As shown in Ref. [16], solving the Dyson-Schwinger-Gorkov (DSG) 



equation and using the instanton packing fraction N/V ~ (200 MeV)^, one can compute 
Mq and A as functions of yU by minimizing the partition function with respect to A. In the 
present work, we consider only the pure NG phase, in which the magnetic susceptibihty is 
well defined in terms of the chiral condensate as shown in Eq. ([T]). We ignore, for simplicity, 
the metastable mixed phases. As a result, the value of the critical chemical potential (yUc) 
was determined to be about 320 MeV. We note that it is about 5% different from that given 
in Ref. [16] in which the renormalization scale is taken to be slightly different from ours, 
A ^ 1/p ^ 0.6 GeV. 



III. QCD MAGNETIC SUSCEPTIBILITY AT FINITE 



In this Section, we compute the VEV of the tensor current of Eq. ([T]) in the the background 
EM field in the presence of the finite /i and provide numerical results for the magnetic 
susceptibility x- order to deal with the background EM field, we employ the Schwinger 
method [6l|29]. Then, we can write the covariant effective partition function with an external 
tensor source T as follows: 

2 



dX Dip Dip^ exp 



■ (10) 



where we have written only the relevant terms to compute the VEV. iD^ indicates the 
covariant derivative, id^ + CqV^, in which Cq and are the quark electric charge and external 
photon field, respectively. The a^i, is the well-known spin antisymmetric tensor defined as 
^tiv = iin^ilv — Ivli^l'^-, and o ■ T denotes a^^T^^" . Differentiating Eq. (10) with respect to 
the external source T, we can evaluate the VEV as follows: 



1 91nZ,ff[T,/i] 



VN( 



dT, 



/d'^p 
-— Trc,T,[^(P,/i)a^^], 



T=0 



{2n 



(11) 



where P stands for p — CqV. Equation (11) can be evaluated further straightforwardly by 



expanding with respect to the electric charge to order 0{eq) as follows: 



SiP, fx) a 



p + iM{P) 



a 



"^■^"''^•'"^%-.M,p)K, + 0«). (12) 



[p2 + M2(p)]2 



In evaluating Eq. (12), we have used the simplified notation p = p + ifi and the identity 
[P^,Pi,] = iCqF^y. The commutation relation in Eq. (12) reads: 



[/',M(P)] = 74P„M(P)] = -2ze,M'(p)(p^F^, + z/iF4.)7.- 



(13) 



Note that there appears a pure electric part proportional to F^^, in the above equation, due 
to the breakdown of Lorentz invariance at finite fi in addition to the usual EM field-strength 
tensor, F^,^. The M' designates the derivative of the effective-quark mass with respect to 
the momentum squared: 



M'ip) 



dM{p) 

Qp2 



(14) 



For convenience, we call the terms proportional to M' as a nonlocal contribution, since they 
arise from the nonlocal (derivative) interaction, whereas the terms without it a local one. 



Considering all the ingredients discussed so far, and performing the trace of Eq. (11) over 
the colour and Dirac spin spaces, we obtain as follows: 



Tr^^fSa 



AN,e,M{p) 
[p2 + M2(p)]2"'''^ 

AN,e,M{p) 
[p2 + M2(p)] 



F 



8N,egM'{p) 

[p2 + M2(p)]2 

8N,egM'{p) 

[p2 + M2(p)]2 



[p^^PaF„u - (/^ ^ l^)] 



magnetic+electric 



SN.egM'ip) 

[p2 + M2(p)]2 



fi^{S^4F, 



4(1 



(15) 



electric 



Here, we have ignored the terms proportional to 0{p) in evaluation of Eq. (15), since they 
become negligible according to the integral identity / d'^pp^ifip^) = 0. 



First, we turn off the pure electric part in Eq. (15) for simplicity, resulting in the magnetic 



contribution to the magnetic susceptibility x, assigned as xm- Incorporating Eqs. (11) and 



(15), and then equating with Eq. ([T]), we finally arrive at the following compact relation for 



dp 
(270^ 



[p2 + M2 



(16) 



One can easily see that Eq. (16) is a generalized expression for Eq. (52) given in Ref. [B] 
by replacing p ^ p + ifi. Switching on the pure electric part, we can write the electric 
contribution to the magnetic susceptibility as xe with an additional nonlocal term, coming 
from the breakdown of Lorentz invariance, as follows: 

d^p \ iJ,^M'{p) 



{2n) 



[p2 + M2(p)]^ 



(17) 



Breakdown of Lorentz invariance 



In Fig. [l| we present the numerical results for the Xm,e(^V'V) (l^ft panel) and the Xm,e 
(right panel) as functions of n, respectively. The magnetic contribution is drawn in the solid 



curve, whereas the electric one in the dashed one. As for the Xm,e(^V'V) depicted in the 
left panel, we show each contribution (local, nonlocal and total) separately. Note that the 
dominant local contributions are the same for the magnetic and electric contributions as 



expected from Eqs. (16) and (17), and starts to decrease beyond /i ~ 200 MeV. At /i = /ic, 
its strength is reduced to about a half of that at yU = 0. On the contrary, the nonlocal 
contributions are almost flat for the NG phase, and show small difference between the mag- 
netic and electric ones, caused by the breakdown of Lorentz invariance, being different from 
the pion weak decay constant [121 [ISl HI]- As a result, the total contributions resemble the 
local one. Moreover, the first-order magnetic phase transition occurs at the critical chemical 
potential, following the chiral restoration, as shown in Fig. [T} Several model calculations, 
using the nonlocal chiral quark model |6j and QCD sum rules [3llll|5], computed x(^^V) cit 
fi = 0. We note that our value ~ 46 MeV is in good agreement with those in Refs. [31 HI [5]. 

Finally, in the right panel of Fig. [T| we show the numerical results for the Xm,e in the 
same manner as in the left panel. Since the chiral condensate is rather flat for the NG phase 
with respect to /i as shown in Refs. [121 [IS], the overall tendency of the Xm,e is very similar 
to that of the xm,e(^V'V)- Again, we observe small difference between the xm and xe, and 
the flrst-order magnetic phase transition at the critical point. Note that in the vicinity of fic 
the Xm,e drecreases by about a factor of two in magnitude compared to its value at /i = 0. 
As mentioned previously, since the Fermi surface begins to be fllled just beyond fi ~ 300 
MeV, which corresponds approximately to the normal nuclear density po ~ 0.17fm~'^, we 
expect that there is a rapid change of the x near po- In other words, the response of the 
QCD vacuum becomes weak but unstable to the external EM source as the density is closed 
to Po, while it remains relatively stable otherwise. 

In Ref. [11], the colour-symmetric one-gluon exchange (OGE) was employed to investigate 
the X in medium in addition to the static screening effect, and was found that a sharp 
divergence takes place in the x near the normal nuclear density, although it depends on the 
screening and the current-quark mass. This observation may be consistent with the present 
results showing the flrst-order magnetic phase transition. 

In fact, partial results of the magnetic susceptibility at flnite quark-chemical potential has 
been already studies within the same framework \l5l. However, there are several important 
points that were not considered in Ref. [15]. Firstly, the breaking of Lorentz invariance 
was assumed to be small, so that it was ignored. Secondly, the chemical potential was 
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H [Me VI /i [Me VI 

FIG. 1: XM,E(i^^V') (left panel) and XM,E (right panel) as functions of /x. The magnetic contribution 
is presented in the solid curve, whereas the electric one in the dashed one. The vertical dashed 
lines indicate the /^c ~ 320 MeV. 

treated perturbatively, so that the calculation was performed by expanding it for numerical 
simplicity. Unfortunately, this naive expansion has led to the wrong behavior of the magnetic 
susceptibility, in particular, as /i is getting close to its critical value. Finally, the finite value 
of the current quark mass was taken into account in Ref. [I5]. However, we soon realized 
that a mere extension of the model by Carter and Diakonov [I6] seems to be unjustified. 
Because of this reason, we are restricted ourselves to consider the chiral limit in the present 
work. 



IV. SUMMARY AND CONCLUSIONS 

In the present work, we have investigated the QCD magnetic susceptibility x induced 
by the external constant electromagnetic field at the finite quark-chemical potential, i.e. 
/i 7^ at T = 0. We focussed on the Nambu-Goldstone (NG) phase, in which the x is 
well defined. We employed the /i-modified nonlocal chiral quark model derived from the 
instanton vacuum at the normalization scale A = p^^ 600 MeV, corresponding to the 
phenomeno logical instanton parameters ^ ^ 1 fm and p ~ 1/3 fm. 

We started with the //-modified Dirac equation to obtain the quark propagator in the 
instanton ensemble. The /i-dependent effective quark mass can be obtained from the Fourier 
transform of the quark-zero mode solution. However, we employed the parameterized dipole- 
type form factor instead of it to circumvent numerical difficulties. The effective partition 
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function was constructed in such way to reproduce the quark propagator in the presence 
of the instanton background. Solving the Dyson-Schwinger-Gorkov equation, as done in 
Ref . [15] , we were able to explore the phase structure through two order parameters, Mq and 
A, corresponding to the Nambu-Goldstone and colour-superconducting phases, respectively. 
As a result, the critical quark-chemical potential /ic was determined to be about 320 MeV. 

Using the low-energy effective partition function with the finite quark-chemical potential, 
we computed the magnetic susceptibility of the QCD vacuum x- We found that the x 
remains stable up to /i ~ 200 MeV, and drops then drastically. At the critical quark- 
chemical potential, fic ~ 320 MeV, the strength of the x is decreased by about a factor 
of two, compared to that at = 0, and the first-order magnetic phase transition takes 
place, corresponding to the chiral restoration. From these observations, we conclude that 
the response of the QCD vacuum becomes weak (insensitive) and unstable to the external 
electromagnetic source near the normal nuclear density in comparison to that for the vacuum. 
Moreover, the effect of the breakdown of Lorentz invariance, caused by the finite fi, on the 
X turns out to be small, Xm ~ Xe'- We note that this tendency is rather different from 
the pion weak decay constant. Related works including explicit flavour SU(3) symmetry 
breaking effects as well as meson-loop corrections at finite quark-chemical potential are 
under progress. 

We want at this point to mention that we have calculated the magnetic susceptibility in 
the NO phase only, that is, below the critical quark-chemical potential. Note that there is, 
however, one caveat: In principle, we can compute the VEV in Eq. ([T| beyond the critical 
quark-chemical potential, i.e., in the CSC phase. We have to keep in mind that the SU(3) 
colour symmetry is broken down to SU(2) x U(l) in the CSC pahse, which causes a lift 
of colour degeneracy. Thus, the quark propagator must be separated into two different 
quark propagators of which each consists of 4 x 4 propagators in the chiral L-R basis [T6] . 
The form of the propagators have been obtained by solving the Schwinger-Dyson-Gorkov 
equation with a systematic expansion in parameters and p/R. Taking consideration 
of these propagators, we can compute the magnetic susceptibility beyond the NG phase, 
though the calculation is quite involved. 

Last but not least, the magnetic susceptibility defined in Eq.Q is not the magnetization 
in medium as done in Ref. [TT]. The magnetization can be rewritten as the VEV of the 
following operator {iq^'^^O'^yq) p- This is another interesting quantity which we can compute 
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within the present scheme. The related investigations are under way. 
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